Restoration of Symmetry by Interactions and Nonreliability of the Perturbative 

Renormalization Group Approach 
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We discuss examples of (l+l)-dimensional models where the perturbative renormalization group 
(RG) indicates a tendency to restore the symmetry in the strong coupling limit. We show that such 
restoration does occur sometimes, but the perturbative RG cannot be reliably used to detect it. 
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Enlarged symmetry phenomenon at a fixed point corre- 
sponds to a situation where the microscopic Hamiltonian 
has a lower symmetry than the fixed point Hamiltonian. 
Any deviations from the full symmetry behaviour scale 
to zero under the renormalization group (RG) flow when 
the critical point is reached. Such phenomenon occurs, 
for instance, in frustrated spin systems near two dimen- 
sions. In the continuum limit, the microscopic Hamilto- 
nian associated with Heisenberg spins on the triangular 
lattice is a non linear sigma model with a 0(3)(Xi 0(2) 
symmetry A RG study of this model in d = 2 -I- e 
dimensions show that the symmetry is dynamically en- 
larged at the stable fixed point to 0(3)(g) 0(3) ~ 0(4) 
The phase transition of canted Heisenberg spin sys- 
tems belongs thus, near two dimensions, to the = 4 
Wilson-Fisher universality class. 

A similar enlarged symmetry scenario has been pro- 
posed by Zhang in the context of high-Tc cuprates 
The author suggests that the phase diagram of these com- 
pounds can be deduced from an SO (5) symmetry that 
unifies antiferromagnetism and D-wave superconductiv- 
ity. It was also argued that even though the SO (5) sym- 
metry is only approximate for microscopic models such 
as the Hubbard or t-J models, it becomes exact under 
the RG flow toward a bicritical point. Microscopic elec- 
tron models with an exact SO (5) symmetry has been 
constructed |^ and this symmetry occurs in the low en- 
ergy sector of two-chain Hubbard systems . In par- 
ticular, Lin, Balents, and Fisher ||6[, using a one-loop 
perturbative RG approach, predict that at half filling 
the weakly-interacting two-leg ladder has an exact SO (8) 
symmetry at low energy and belongs to a D-Mott phase 
where the short range pairing correlations have an ap- 
proximate D-wave symmetry. The crucial point of their 
analysis stems from the fact that some coupling constants 
of the microscopic Hamiltonian flow to strong coupling 
but with ratios converging to fixed values and the effec- 
tive model corresponds to the SO (8) Gross-Neveu model 
in the strong coupling limit. 

The question that we shall adress, in this letter, is 



whether one can reliably deduce restoration of a symme- 
try by interactions in the strong coupling limit from weak 
coupling RG computations. In the following, we shall 
consider several (l-l-l)-dimensional models using a non- 
pertubative approach to investigate the possible dynami- 
cally enlarged symmetry in the strong coupling limit. Let 
us start with the U(l)-symmetric Thirring model defined 
by the following Hamiltonian density: 



a=l 



(1) 



with gi = g2 — g±, and g^ = g||. The (iso)spin cur- 
rents are expressed in terms of the fermion fields (^a) 
as: J° = (l/2)^Q7^cr°^^/3, being the Pauh matrices 
and our Euclidean conventions for the gamma matrices 
are: 7° = cr^, 7^ = cr^. The Thirring model (Q) is exactly 
solvable by the Bethe ansatz and the solution was care- 
fully studied . This gives us the first opportunity to 
test whether restoration of symmetry occurs and if the 
one-loop approximation in the RG equations can serve 
as a reliable tool to detect it. The RG equations for this 
model are of the Kosterlitz-Thouless form: 
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where the dots indicate derivatives with respect to the 
logarithm of the length scale: g = dg/d\nl. The phase 
portrait derived from these equations is presented on Fig. 
1; there are regions where the RG trajectories flowing 
to strong coupling converge giving appearence that the 
original U(l) symmetry is enlarged to SU(2). The strong 
coupling regime can be separated into two regions. The 
first one corresponds to g^^ > g± > (marked A on 
Fig. 1) where the Thirring model is equivalent to the 
sine-Gordon model. We shall be interested only in the 
region where the symmetry is apparently restored. In 
this region the spectrum of the sine-Gordon model con- 
tains only kinks and anti-kinks with equal masses and no 
bound states. All information about interaction between 
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these particles is contained in their two-body S-matrix 
which depends explicitly on the renormalized anisotropy: 

■■ — arccos[cos (7||/ cosg^J f 



(3) 



where the latter equality holds at small bare couplings. 
The S-matrix is an analytic function of ^ at ^ — > and 
corrections to the isotropic limit are order of ^ at small 
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FIG. 1. Flow diagram of the U(l) Thirring model. 



It turns out, however, that restoration of the symmetry 
is much stronger in the other region, namely g± > \g\\\. 
Following Japaridze, Nersesyan and Wiegmann (JNW) 
who studied it we shall call it region C (see Fig. 1). 
The JNW solution can be generalized for higher repre- 
sentations of the SU(2) group and also for other semisim- 
ple groups All these integrable models are defor- 
mations of the Wess-Zumino-Novikov-Witten model by 
marginally relevant perturbations. Let operators J" and 
J° be the left- and right currents of a semisimple group 
Q satisfying the level k Kac-Moody algebra: 

[r{x), J\y)] = ^S'^'S'ix -y) + ir'^,r{y)5{x - y) 

(4) 

with a similar relation for the right spin current (J) 
whereas J and J commute each other. In Eq. (Q), /'^^'^ 

for 
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are the structure constant of the algebra (/"' 
the SU(2) group). The corresponding Hamiltonian den- 
sity can be written in terms of the currents (Sugawara 
form) : 

2„ G _ _ G 

n = - — ^[: rr ■. + : rr ■] + ^ gabJV' (s) 



a=l 



where G is the number of generators of the group Q and 
the Casimir operator (c„) in the adjoint representation 
is defined by: c^Saa = f'^^'^f'^^''. For each Q and k there 
is at least one single-parametric family of coupling con- 
stants (we do not consider the common factor as a 



parameter) at which the model is integrable. In particu- 
lar, for 0(2n) groups there are two families of solutions. 
There is one trouble, however, - for all groups higher than 
SU(2) the parity is broken: gab gba (to get an idea one 
may look at Refs. |pO|-p^ where some related Hamiltoni- 
ans are written exphcitly). For the SU(2) group there is 
no trouble and the model ^ becomes a generalization of 
the Thirring model with two coupling constants g± and 
g|| . At fc = 2 the isotropic SU(2)-invariant version of the 
model (||) becomes the 0(3) Gross- Neveu model. This 
model with = is equivalent to the model suggested 
by Egger and Gogolin in the context of theory of 
carbon nanotubes: 

i ^ 

'H^ ^{rad^ra - Ldxla) + g r3h{rili + r2^2) (6) 

a=l 

where ra,la are right and left moving real (Majorana) 
fermions. The equivalence between the SU(2)-invariant 
version model ^ with k = 2,g^\ =0 and model (||) 
stems from the fact that the SU(2)2 spin currents can 
be written in terms of three real fermions [p^ljl^ ] : J° = 
— ie'''"^ZfcZc/2 with a similar relation for the right-moving 
current. The RG equations of the model (^) in the first 
loop approximation does not depend on k. Thus in the 
region of small bare couplings the phase portrait of the 
model is still the one presented on Fig. 1. 

The Bethe ansatz equations for the SU(2)-based model 
(^ were obtained by Wiegmann and Kirillov and 
Reshetikhin [0 and are given by 

M 

[ekixa + f/2fi)ek{xa - f/2^i)f = J] 62(0;, - Xb) (7) 

6=1 

where 2N is the number of fermions, M is the number of 
up spins. In the region A the functions e„(a;) are hyper- 
bolic whereas in the region C they are trigonometric: 



sm[ii{x — i«-/2)] 
sin[/i(a; + in/2)] 



In the region C the parameters /i and < / < tt are re- 
lated to the bare coupling constants (which are supposed 
to be small) g: 



-gl + gl, cot / = g\\/^-gj + gl. 



(8) 



The results of the previous analysis show that in the re- 
gion A the symmetry is not restored in the model (|^). 
The region C for the case fc > 1 has never been studied, 
but it is easy to see that the results obtained in Ref. j|] 
for the case k = 1 are straightforwardly generalizable for 
greater fc's. In the region C the models (0j|) have the 
spectral gap m « Aexpf— 7r(7r — /)/2/i]. The important 
result obtained in Ref. g is that in the strong coupling 
limit the anisotropy of the original Hamiltonian leads 
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only to exponentially small corrections to physical quan- 
tities ~ exp(— 7r^//Lt). From the technical point of view 
the most important difference between the Bethe ansatz 
equations in the two regions comes from the fact that in 
the sine-Gordon region they contain hyperbolic and in 
the region C - trigonometric functions. In the latter case 
the Thermodynamic Bethe Ansatz (TBA) equations dif- 
fer from the equations at the isotropic point (/i — *■ 0) only 
by the fact that real parts of the spectral parameters Xa 
belong to a finite interval < 3iex < 2tt/ n and the ker- 
nels in the integral equations are periodic with the period 
27r//i. As a matter of fact, one can obtain TBA for the 
sector C from the isotropic TBA replacing the kernels by 
their periodic generalizations: 



Kc{x) ^ ^ Kisotrix - 2Trl/^i). 



(9) 



At temperatures comparable with the spectral gap the 
thermodynamics is dominated by a small part of the 
available parameter interval and the difference between 
periodic and isotropic kernels becomes insignificant (ex- 
ponentially small in n/fi^ to be precise, because kernels 
in TBA equations are exponential functions of x). 

As we have mentioned, one application of this result is 
the Egger-Gogolin model of the carbon nanotubes. In the 
strong coupling limit, one can conclude that the model 
is equivalent to the 0(3) Gross-Neveu model, which 
can be described using the results obtained by one of the 
authors There is another application: the model 

(^ with k ^ oo corresponds to the model of classical 
Heisenberg antiferromagnet on a triangular lattice. At 
T = 0, as consequence of frustration, the system orders 
in a spiral state with the order parameter consisting of 
two mutually perpendicular unit vectors. In the contin- 
uum limit, the long distance properties are described by a 
nonlinear sigma model with an order parameter in SO (3) 
[0,^. Using a SU(2) matrix representation, this model is 
described by the anisotropic Principal Chiral (PC) model 
with the energy density: 



£ = 



2 



(10) 



where the bare coupling constants verify 2p±^ = pn > 0. 
The vielbein 17° is expressed in terms of the matrix g 
belonging to the SU(2) group by: fl'^ — ilr{a°'g~^d^g). 
Since p±_ ^ p\\, the bare action has a SU(2)xU(l) sym- 
metry. Using the Polyakov- Wiegmann approach |19| , one 
can relate the model (10) to a generalized Thirring model 
considered above and thus study the possible dynami- 
cally enlarged symmetry of the model ( |lO| ) in the strong 
coupling limit. This correspondence can be established 
by considering the partition function of the anisotropic 
PC model: 



Z = jvnid (F^,) cxp (^-1 1 £x Pa^l^l^ (11 



where the delta constraint indicates that the gauge field 
crT2^ = —2ig~^d^g has a zero field strength. One has a 
fermionic representation of this constraint with the iden- 
tity |9|: 

5{F^,) = lim / V^.^V^,^ 

exp J d^xJ2^,„ {-/""S^pd^ + i7^f^>S/3) ■ (12) 

Inserting this result in (^ij) and integrating over the fi" 
fields, one obtains the Hamiltonian density of a general- 
ized Thirring model: 



fc 3 



(13) 
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where the color index k tends to infinity and the current 
is given by: = i'^/'2)T,^=l^^c.l^<y^p'^^p■ Since the 
bare coupling constants of model (|l3|) verify pj^ > p'^^, 
the anisotropic PC model (|l0| ) belongs to the region C of 
the phase portrait of the SU(2)-invariant model (||) with 
k — > oo. Therefore, according to the arguments given 
above, the anisotropic PC model acquires an enlarged 
SU(2)xSU(2) symmetry in the strong coupling limit. In 
this limit its solution coincides with the solution of the 
isotropic Principal Chiral Model (see Refs. [p^ , pO| ). 

The third example we consider is directly related to the 
work of Lin, Balents, and Fisher It is the anisotropic 
Gross-Neveu model with the Hamiltonian density {U > 
0) 

JV 
a=l 

-giii^oi'of + giMo){'4>a^a) + Ui^Pa-^Paf ■ (14) 

One can obtain the RG equations of the model (jlj) : 

^ ^ _2_iN-l) ^ l_ 
27r 
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(15) 



These equations indicate a restoration of the symmetry 
in the strong coupling limit: the model has an isotropic 
behaviour at strong coupling with U changing its sign and 
flowing to the strong coupling. We shall see however, that 
this isotropy is not reflected in the excitation spectrum. 

One can, first, bosonizc the Dirac fcrmion iJjq with the 
introduction of a bosonic field <& to get another formula- 
tion of (0): 



N 

E 

a=l 



i^ali^id^^lpa - gC0S/3$(V'aV'a) 

+ U{Ma? (16) 
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where g '-^ g and for gi > the scahng dimension of the 
cosine-operator is c? = 2/?^ < 1. The operator containing 
the cosine term is then relevant and one may think that it 
wiU generate a mass gap for aU branches of the spectrum. 
Thus the mass gaps in this case would be generated by 
the repulsive interaction in the 0-channel; this mechanism 
was discussed by Shelton and one of the authors in con- 
text of doped ladder models in Rcf. where a model 
very similar to ( p^ ) with N ^ 2 was considered. The 
model ( p^ has 0(2)xO(2N) symmetry. According to 
the one-loop RG equations, this symmetry is enlarged at 
strong couphng to 0(2N -|- 2). Let us consider the model 
( |l6| ) in the large TV limit to test this possible restora- 
tion of symmetry. In this case we can use the mean-field 
approach to decouple the interaction between and a- 
sectors: H ~ Hsg + Ti-oN with 

HsG = 7^(9^*)' - Mcos/3$ (17) 

N 

'HgN = ^ i>aJtidfi1pa - moii'ai'a) + U (ipall^a)^ (18) 
a=l 

where M — g < {■ipa'4'a) > and mo ~ g < cos/3$ >. 
Using the exact expression of vacuum expectation val- 
ues of exponential fields in the sine-Gordon model [p2[ , 
one can relate rriQ to the soliton mass mgQ of the sine- 
Gordon model ([l7|): TOq = gCd rrigQ where the numerical 
prefactor Cd is known [p2|: 



nT{l-d/2) 



2(2 - d) sin[7rd/(2 - d)]T{d/2) 



[TTd/2{2 - d)] 



2^ 



r[l/(2-d)]r[l-d/2(2-(i)] 



Using the large N limit of model ([Tsl), one obtains the 
following expressions: 



M 



~gN 



2tt 
Too 



TOln(A/TO) 



1 -f- [2UN/Tr] ln(A/TO) 



(19) 
(20) 



In order to close the different equations of the large N 
limit of the model (^6|), we need a relation between M 
and msG which is given by: M ~ DdmtTf with p3| 



2r(d/2) f V^T[l/{2-d)] \ 
7rr(l - d/2) \ 2r[d/2(2 - d)] J 



2-d 



We shall assume that in the limit N ^ oo, UN ^ 
U*,g'^N'^ = g* remain finite and that U* ln(A/TO) >> 1. 
In that case, we obtain the mass spectrum of the model 
in the large N limit: 



msG 



TO ln(A/TO) 



[aU* Dd) 

2^Z?45*)i/(2-2'i)(c<i/4Z?dC/*)(2-'i)/(2-2<i). 



As wee see, the ratio of the masses msG/'m ~ VN is 
large. Therefore, this result indicates that the symme- 
try is not restored by interactions in the strong coupling 
limit. 

From our discussion of generalized Thirring models, we 
conclude that there are two possible sectors where the 
symmetry is apparently restored in the strong coupling 
regime. In the first one, called the sine-Gordon sector, 
the parameter characterizing anisotropy does not renor- 
malize to zero at strong coupling. Its value depends on 
the bare coupling constants and is unlikely to be small 
outside the region of weak bare interactions. There is 
another region, however, the region C, or the trigono- 
metric region where the influence of the anisotropy is 
exponentially small. Here the restored symmetry has 
more chances to survive outside of regions of weak bare 
couplings. Both regions arc indistinguishable within the 
perturbative RG approach in the infrared limit. We 
have also given an example where the RG equations in- 
dicate isotropy, but the nonperturbative approach based 
on large N-approximation shows stricking anisotropy. 

We are grateful to F. Essler, A. Gogolin, A. Nersesyan, 
and S. Sachdev for interesting discussions and interest to 
the work. 
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